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Lorentz violation caused by the spacetime defects, which are expected to reveal themselves in 
the photon sector. A modification of the photon dispersion relation may result as a quantum 
effect and we compute it at leading order perturbation theory. 

The outcome of the calculation is that the photon dispersion law remains conventional, 
if the defect distribution is dense, homogeneous, and isotropic. This result sheds some new 
light on Lorentz violation in the framework of a small-scale structure of spacetime. It shows 
that Lorentz invariance can be preserved even in the presence of a spacetime structure that is 
supposed to emerge at the Planck scale. This conclusion has already been drawn on general 
grounds in other publications, where the current paper delivers a demonstration by a direct 
computation in a simple model. 
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I. INTRODUCTION 

Physics at the Planck scale is considered to be the holy grail in present fundamental research. 
Up to now it has not been rigorously shown which physical phenomena occur at the Planck scale 
and how they can be described within a mathematical formalism. For this reason it makes sense 
to construct simple models in order to mimic effects expected to occur at this scale. The property 
"simple" means using underlying concepts of established theories such as general relativity and 
quantum theory, which are well understood and hold for energies much smaller than the Planck 
energy. 

One fundamental concept are Einstein's field equations in general relativity linking energy 
density to spacetime geometry. Another one is Heisenberg's uncertainty principle in quantum 
theory saying that two complementary particle properties are always endowed with uncertainties. 
Assuming that both still hold at the Planck scale, energy uncertainty will result in an uncertainty of 
spacetime geometry or even topology. Thus, under this assumption, spacetime metric coefficients 
and spacetime curvature will begin to fluctuate at the Planck energy. Such fluctuations are often 
referred to as spacetime foam or spacetime defects [THS]. 

There exist various models for spacetime defects with nontrivial topology. One way in which 
to obtain a spacetime defect is to cut an open set out of Minkowski space and to impose certain 
conditions on the remaining boundary. This then leads to a spacetime Ai having a different 
topology. For example, cutting out an open ball of M 3 and identifying antipodal points on the 
boundary results in M. = IR x (IRP 3 — {point}), where IRP 3 is the three-dimensional orientable real 
projective space [7|. 

However, the interest not only lies in the defects themselves, but in their influence upon the 
propagation properties of particles, for example photons. In jS], the modification of the photon 
dispersion law is investigated for certain classes of defects. The method used is to consider the 
scattering of an electromagnetic wave at one single defect. Certain conditions for the physical fields 
are then set on the boundary of the defect and Maxwell's equations are solved by introducing a 
correction field. Even for a single defect this is a difficult task and the approach would be more 
difficult for two defects and impractical for a large number of defects nearby. 

Since we are interested in the propagation of photons through a spacetime foam made up of many 
defects, we proceed with an alternative possibility that was initiated in [9j. Here the CPP-anomaly 
is derived for a non-Abelian gauge group 5*0(10) on two spacetime manifolds with nontrivial 
topology: a spacetime with a linear defect M = I x (R x (IR 2 \ {0})) = IR 4 \ IR 2 and a spacetime 
with two identical open balls removed from IR 3 and points on their boundary properly identified 
(wormhole). It is shown that the CPP-anomaly occurs for the Abelian subgroup U(l) C SO(10), 
as well. The anomaly gives a contribution to the effective action as an PP-term with the field 
strength tensor P and its dual P. Since it is tremendously difficult to obtain the anomaly for 
several defects, the introduction of a background field is mandatory. This field does not describe 
any microscopic defect properties and, hence, serves as an effective approach for the case when the 
photon wave length is much larger than the defect size. 

We will follow this idea and describe a single defect as pointlike, where it is assumed to be 
time-dependent — contrary to Ref. [9]. Such defects are distributed randomly in Minkowski space 
resulting in an effective "random" background field. Furthermore, the distribution of defects is 
taken as being Lorentz- invariant. We study whether and how the dispersion relation of photons is 
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affected by such a time-dependent and Lorentz- invariant background. 

The outline of this paper is as follows. In Sec. [TT] we introduce the action of the effective theory 
that forms the basis of the article. This action describes the interaction between photons and 
the defects that is mediated via a real scalar field. Section |III| gives a description of how to treat 
photon propagation through a distribution of many pointlike defects that are put randomly at 
distinct points in Minkowski spacetime. In Sec. |IV| the focus will be on the perturbative solution 
of the photon field equation. In the first part the solution is obtained by inserting a perturbative 
Ansatz into the field equation. In the second part we demonstrate how the same result follows from 
a diagrammatic approach. As a next step we compute the leading-order solution of the photon 
field equation in Sec. |VJ where in this context a renormalization procedure has to be performed. 
In the follow-up section, VI the scalar field equation is solved to leading order in perturbation 
theory as well. Both results are combined to calculate the modified photon dispersion relation in 
Sec. VII and the physical meaning of the result is then discussed. Throughout the paper certain 
assumptions will be made so that the calculation is feasible. In Sec. VIII we make a couple of 
remarks on how the result may change if certain assumptions are dropped. In the penultimate 
section, |IX[ we go on a brief excursion to /T-symmetric quantum field theory in the context of 
the special model proposed. The last section, |Xj gives a summary and a conclusion on the results. 
The most important computational steps can be found in Apps. [XJ [BJ [Cj and [Dj 

Throughout the paper natural units are used with h = c = 1. For some occurrences, c and h 
will be reinstated for clarity. 



II. ACTION OF THE EFFECTIVE THEORY 



We wish to describe photon propagation through a Lorentz-invariant distribution of time- 
dependent, pointlike spacetime defects. The action upon which our considerations are based is 
given as follows [10J: 
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where F fMU (x) = d^A u (x) — d u An(x) is the field strength tensor of the U (1) gauge field A^{x). The 
first term in the action is the standard kinetic term of the free photon field and the second fixes 
the gauge (we use Feynman gauge). The third contribution contains the kinetic and mass terms 
of a free real scalar field 4>{x), where is a parameter with mass dimension —1. The scalar field 
<p(x) itself has mass dimension zero and its bare mass is given by l/i>'°). 

The fourth expression involves the pointlike spacetime defects sitting at distinct spacetime 
points x n . The defects are effectively described as randomly distributed single four-dimensional 5- 
functions, each carrying some charge e G { — 1,1}. Concretely, this means that each defect appears 
at a single point in three-space for an infinitesimally short amount of time before disappearing 
again. This illustrative picture corresponds to what a theorist has in mind when thinking on a 
simple spacetime foam. 
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TABLE I. Fields and nondynamical objects (defects) appearing in the action of Eq. (2.1 1 with corresponding 
mass dimension, bare coupling constant A' ) <C 1 and e = ±1. 



Only the field (f)(x) is assumed to couple directly to the defects via the charge e. Finally the 
interaction of the photon field with the defects is described by the last term, where the interaction 
is mediated by the scalar field eft. The latter term involves the dual field strength tensor 



F^(x) = \e^F g(T (x) 



(2.2) 



with the four-dimensional Levi-Civita tensor e> iv s a ' . The motivation behind such an interaction is 
that it appears in the effective action in the context of the CPT-anomaly [9]. The function f((ft(x)) 



in the last term of Eq. (2.1 ) can be arbitrary, in principle, but it is assumed to be sufficiently well- 
behaved. In order to keep the model simple we choose f{(j){x)) = (p(x). Note that the fifth and 
sixth terms explicitly break gauge invarianceQ All fields are defined in Minkowski spacetime with 
metric r]^ = diag(l, — 1, — 1, — 1). To summarize, the field content of the theory is presented in 
Table HI 

If we wanted to couple photons to a conservecj^] fermionic current j^(x) = ip(x) r y' J "ifj(x) with 
the Dirac field ip{x) and standard Dirac matrices 7^, the modified theory could be coupled to the 
Dirac theory of standard spin-1/2 fermions of charge e and mass M: 



d A x ^(x){ 7 ^[ia M - eA^x)] - M}V(3 



This would make the complete action of the theory 

irac > 



(2.3) 



(2.4) 



with S e s given by Eq. (2.1) and Sbirac by Eq. (2.3). The description of the spacetime foam 
model by the action given is the fundamental assumption of this paper. It will be referred to as 
Assumption (1)0 



III. STATISTICAL TREATMENT OF A LARGE NUMBER OF SPACETIME DEFECTS 
A. Distribution of defects in Minkowski spacetime (sprinkling) 



We intend to distribute spacetime defects in four-dimensional Minkowski spacetime Ai in a 
Lorentz-invariant manner. This will be possible if defects are distributed according to a "Poisson 

1 They are not invariant under a gauge transformation of the field 4>, namely 4>{x) i-> — <j)(x). To be crystal clear, 
this gauge transformation has nothing to do with the U(l) gauge transformation of the photon field. 



2 By quantum corrections the explicit violation of gauge invariance in the action S c ff of Eq. (2.1 1 may give rise to 
an anomalous nonconserved current. However, this effect is expected to be suppressed by (A^) 2 . 

3 In what follows, several further assumptions will be taken. In such a context the word "assumption" is abbreviated 
as "Ass." in combination with a number and optional small Latin letters. 
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process" (i.e. a sprinkling). The result of the Poisson process is a Poisson distribution of defects 
throughout the spacetime. This means that the probability of observing n defects in a rectangular 
spacetime region with side length ££ and volume 



V = / d 4 x J-det^) = / d 4 x = ^ 4 , (3.1) 



region \x»\<M/2 

^£{0,1,2,3} 



is given by: 



F,(V) s g°M . (3. 2) 



nl 



Herein, g is at first a parameter that characterizes the distribution. Note that Eq. (3.2) is a valid 
description for a probability distribution, which is both isotropic and homogeneous. The most 
natural choice for a distribution of pointlike — i.e. zero-dimensional — spacetime defects is an 
isotropic one if we do not take into account any mechanism producing defects that make space 
anisotropic (for example, defects similar to cosmic strings [II]). When the volume of the region 
approaches zero (V i— > 5V with an infinitesimal value SV) , the probability of finding a single defect 
in that region is proportional to the volume 

P n=1 (5V) = g5V + 0{5V 2 ) . (3.3) 
On the other hand, the probability of finding more than one defect is negligible: 

P n>1 {5V) = 0(5V n ). (3.4) 
An explicit realization of the Poisson process is then described by the following steps [12]: 

1) Divide V into small boxes with volume V. 

2) Then place a defect into each box with probability P = gV . 

3) The Poisson process is obtained in the limit V i— > 0. 

The Poisson process or "sprinkling" is invariant under any volume-preserving linear transformation 
and in particular it is invariant under Lorentz transformations (this happens because the process 
only depends on the spacetime volume V). Moreover, it has been shown in |13] that the realizations 
of the Poisson process are Lorentz-invariant individually as well. Lorentz invariance in this context 
has the following meaning 

"The discrete set of sprinkled points must not, in and of itself, serve to pick out a preferred 

reference frame." |12j 

That is, the statistical properties of the distribution of defects (e.g. the mean density of defects) 
do not depend on which reference frame we choose to measure them in (see Fig. [T]) . 

We emphasize that the number of defects contained in different regions of equal volume V is 
not constant but fluctuates from region to region 



N(V) = (N{V))±SN{V), 



(3.5) 
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(a) = (b) = 0.7 



FIG. 1. Example of sprinkling in a finite region of a two-dimensional spacetime (ID-space x time) as it 
looks in two different inertial frames. The second frame is boosted along the positive spatial axis with a 
Lorentz boost factor /3 = 0.7 with respect to the first frame. This boost changes the shape of the spacetime 
region. However, the mean density (/9 bs) of defects is the same in both frames. As an illustration compare 
the two enlarged regions of both distributions. Note that in the boosted distribution not all defects are 
shown. 



where the mean number of defects and the standard deviation are respectively 



(N(V)) = nP n (V) = gV , 6N(V) 



n=0 



Y j {n-(n)yP n (V) = ^g~V 



(3.6) 



n=0 



These results can be explicitly obtained by using Eq. (3.2 ). Therefore we can identify the parameter 
g with the mean density of defects: {p Q bs) = Q- Moreover fluctuations in the number of defects also 
imply fluctuations in the density /? bs from a region to another 



Pobs = (Pobs) ± Sp 



(N(V)) ± 5N(V) 



V 



V 



Q± 



(3.7) 



Nevertheless, these fluctuations become negligible when the mean volume occupied by the defects, 
namely Vd = is much smaller than the volume V of the region considered 



1 



<C V 



Sp 
Q 



1 



0. 



(3.8) 



Thus we can regard the density as being constant so long as we consider scales that are much larger 
than the mean separation between defects: N{V) = qV = const, when V S> Vd- Being able to 
perform computations with a globally constant volume means that the distribution is homogeneous 
as well. Otherwise the density g would only be defined locally. We will refer to the isotropic and 
homogeneous distribution of spacetime defects as Ass. (2a) and Ass. (2b), respectively. 



B. Derivation of auxiliary functions 



In this section we intend to derive a set of functions that will be extensively needed for the 
solution of the scalar and photon field equations. First of all, we investigate the problem of photon 
propagation through a spacetime with N defects in the finite subset with volume V given by 
Eq. (3.1). The corresponding functions and fields are denoted with indices W or l M\ We now 
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consider the "free" field equation for (f>, i.e. neglecting the coupling to the photon but not to 
the defects. In order to solve the field equation, the truncated Fourier transform of <j>(x) when 
restricted to the box is necessary 

W)= J dWifc.x)«x). (3.9) 

\x»\<&/2 
^£{0,1,2,3} 

Analogously we define the inverse Fourier transform by 

d 4 k 

exp(i&; • x)g A r(k) . (3.10) 



(2vr) 4 

/iS{0,l,2,3} 

In configuration space the coordinates are restricted by the side length of the box considered. This 
corresponds to a minimum value in momentum space, which is manifest in the integration limits. 
In principle, for a finite volume the Fourier transform would correspond to a Fourier series with 
discrete coordinates or momenta. However for simplicity we will assume a continuous spectrum. 
This does not play any role as we will eventually generalize the results to the whole of Minkowski 



spacetime anyway. We use Eq. (3.10) as an Ansatz to obtain the following solution of the field 



equation of 4>{x) in momentum space 

g M {k) = VMH{k)G N {k) , (3.11a) 



H(k) = 



.( 6 (0))S 



k 2 -l/(6(°)) 2 + ie 



G^f{k) = — = ^2 £ i exp(i/c • x { 



(3.11b) 



The solution depends on the bare mass of the scalar field and the distribution of spacetime defects. 
Furthermore it consists of two contributions. The first, Gj^(k), solely describes the defects and the 



second, H(k), serves as a mediator between the defects and the photons. In principle Eq. (3.11a) 
can be regarded as being a solution of the classical field equations restricted to a finite box. Later 
on we will need to consider the product of two such solutions evaluated at different momenta 
- especially for the perturbative photon field. The function g^{k) then serves as an effective 
background that describes the influence of the spacetime defects on the photons. Therefore we will 
refer to it as an effective background field in what follows. 



Since Eq. (3.11a) gives the background field with respect to the number of defects, we define 



the corresponding background field with regard to the volume: 



g a (k) = VVH(k)Ga(k), G M {k) 



^5>exp(i^ 



(3.12) 



i=l 



We keep in mind that 4> couples to the photon via the last term of Eq. (2.1 ). After establishing all 
intermediate results we compute the limits M H > oo and 3$ ^ oo, respectively. 

In the subsequent paragraphs we are going to perform a statistical treatment of photon propa- 
gation through the background field. In light of this we will encounter products of functions defined 
by Eqs. (3.11a), (3.11b) (or Eq. ( 3. 12[ ) ) , where each depends on a different four-momentum. Such 
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products will have to be summed over the total number J\f of defects distributed in the whole of 
Minkowski space. 



According to Eq. (3.2) the defect density does not change with respect to Lorentz transfor- 
mations. However this does not mean that the defect positions are Lorentz-invariant. Therefore, 
besides summing over the number of defects we are also going to integrate over all Lorentz trans- 
formations from zero to an upper limit A'. We introduce this cutoff to avoid difficulties from 
integrating over the noncompact Lorentz group. We are then ready to compute the sum of com- 
plex exponential functions (random phases), each evaluated at the spacetime point of a defect. 
The corresponding result will be needed later: 

A' N A' 

/ dA ^2 expfiA: • (Ax n )] ~ / dA g / d 4 x exp[i/c • (Ax)] 

J n=l ^ ^ 

A' 

= J dAgJ dV | det(A" 1 )| exp[iA; • x'} 

A' 

= Z(2ir) 4 g5 {i) (k), E = JdA. (3.13) 

In the second step we have approximated the sum over n as an integral over x: 

M 

^ ^ / dn= g / d 4 x. (3.14) 

n=l J 

This is possible since the sprinkling procedure ensures the proportionality between the number of 
defects and the volume: dn = gd 4 x. Furthermore the defect distribution is assumed to be dense. 
The latter is an important issue not only for calculational but for physical reasons as well. In Ref. [8] 
a classical spacetime foam with topologically nontrivial defects having a particular size b and a mean 
separation I is considered. Bounds obtained from the absence of vacuum Cherenkov radiation lead 
to the constraint b/l < 10 -7 within the spacetime foam model considered in this reference. Hence, 
for spacetime defects of the size 10 2 x Lp\ with the Planck length Lp\ = Gh/c 3 ~ 1.62 • 10~ 35 m, 
where a classical approach is supposed to be valid, the defects would be separated by at least 
10 9 x Lpi ~ 1.62 • 10~ 26 m. Even if their separation is larger by several orders of magnitude, the 
approximation of a dense distribution still makes sense. We assume that this conclusion can be 
applied to the spacetime foam model with pointlike defects investigated here. This is fortified by 
Table |V| in App.|Aj The dense distribution of defects will be referred to as Ass. (2c). 



The last assumption in Eq. (3.13) when computing the remaining integral is an infinite space- 
time. Current cosmological data implies that we live in a flat universe of finite age. For the 
curvature radius of the universe, which is related to its size, only upper bounds can be given within 
the ACDM model [13] • Therefore it has not as yet been clarified whether the universe has a finite 
or infinite volume. However considering photon propagation on time scales that are much smaller 
than cosmological time scales the Friedmann-Robertson-Walker metric for a flat universe corre- 
sponds to the Minkowski metric to a reasonable approximation. In order to keep the model as 
simple as possible, we do not describe cosmological effects and assume the spacetime volume to be 
infinite. In what follows we will refer to the infinite spacetime volume as Ass. (2d). 

After writing the sum as an infinite spatial integral, we change the variables to x \— > x' = Ax 
and employ the fact that | det(A _1 )| = 1. The integral leads to a four-dimensional 5-function. The 
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defect index 



i 

3 =h i 



sign of charge 



±1 



defect position 



X E ,j 



XE,; 

xea 



Sx 



E.i 



TABLE II. For a large number of defects we consider a single defect at Xe\ with a specific charge. Then 
the probability of finding a neighboring defect at xej with a charge of the opposite sign, where Xej has a 
small Euclidian distance to xe,u is close to 1. 



result does not depend on the total number TV of defects any more but only on the density q. The 
latter acts as a constant of proportionality which is independent of the reference frame. Because of 
this the final result is Lorentz-invariant. The Lorentz transformation of the defect positions does 
not play a role and what remains is the integral over all Lorentz transformations. It corresponds 
to the "volume" 5 of the Lorentz group when restricted to an upper cutoff A' and it is a mere 
number. 

We can now compute the product of two functions Gtf — followed by an integration over all 
Lorentz transformations up to the cutoff A' — in the limit of large M: 



I 



dA lim Gj^(k)Gj^{p) = [ dA lim — 
A/V->oo J AA-»oo J\l 



^exp[i(/c +p)(Axi)} + Pmn 

8=1 m^n 



(3.15) 



The quantity P mn involves the product of charges of different defects: 

Pmn = £ m s n exp[i/c • (Aa? m )] expfip • (Ax n )] . (3.16) 

The sum over P mn has to be evaluated for a large number of defects. This will be done in the 
following few lines. We transform coordinates to Euclidian space via a Wick rotation, where these 
coordinates will be marked by an index i E\ Now consider a small hypercube H an (xE, n ) around a 
defect at xe,u with side length 

y\V 4 f\\ l / A 



a n=[jf) =[-) . n€{l,...,M}. (3.17) 

On average every defect lies within such a hypercube. For N 3> 1, given a defect at x n with charge 
e = ±1, a partner with a charge of opposite sign e = =Fl can be found at a distance 5xe,u that is 
of the order of the side length a n of a small hypercube (see Table |n]) . We can therefore perform 
a Taylor expansion for P mn in the small parameters 5xe,u- For the sum over P mn for m ^ n we 
consider defects at XE, m and xe,u with charges e = 1 and their neighbors at SxE, m and 5xe,u with 
charges e = —1. Substituting x' Ei = Axe,% we obtain 



Pmn = Y £m£n eX P( ik ' X E,m) exp(ip • x' E<n ) 



E 

m^n m^n 



X] K+ 1 ) 2 exp(i/c • a4,m) exp(ip • x' E>n ) 

m^n 

+ (-1) 2 exp[i/c • (x' E m + 5x' E m )} expfip • (x' E>n + 5x' E , n )} 
+ (+ 1 ) • (-1) exp(i/c • x' Ejm ) exp[fp • (x' E , n + Sa/ Efn )] 
+ (-!)• (+1) exp[i/c • (x' EjTn + Sa/ E m )] exp(ip • x' E n )} 
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= 0{5x E ^ m ■ 8x' Ejn ) . (3.18) 

As a result the linear term vanishes and all further contributions are suppressed by small distances. 
Thus for a large number of defects in the sum over all P mn , contributions from neighboring defects 



with opposite charges compensate each other. The second term of Eq. (|3.15|) averages out 
lim —. 

A/Woo J\f 



We obtain the following result for the product of two functions Gj\f, each evaluated at a different 
momentum 

A' 

f dA lim G N (k)G N (p) = Z lim B 5 {i) {k + p) . (3.20) 

J A/V-»oo A/V-»oo J\l 

The integral over all Lorentz transformations up to a boost limit A' can be separated out again. 



Furthermore, the functions G@{k) and gs?,{k) defined by Eq. (3.12) obey a similar relation 



[ dA lim GM(k)G M {p) = E lim i (2tt) V (4) + P) • (3.21) 
J se^roo Vh->oo V 

The physical interpretation of the result obtained is as follows. If the scalar field <j> scatters at 
a defect it can either transfer momentum to the defect or absorb momentum from the defect. 
Averaging over many defects — implying the limits i — ^ oo, V i- > oo or ^ i- >• oo — leads to 
zero average momentum transfer at each defect. This means momentum conservation for the <f>- 
field. Averaging over infinitely many randomly distributed defects results in a translation invariant 
theory that clearly obeys the property of momentum conservation. 

The results obtained depend only on the density of defects and (for some) M or V. Hence these 
can be generalized to the whole of Minkowski spacetime. As mentioned this corresponds to blowing 



up the spacetime region of Eq. (3.1), namely to the limit M \— > oo. 



IV. PERTURB ATIVE SOLUTION OF THE FIELD EQUATIONS 

The dispersion relations of both the scalar field and the photon follow from the appropriate 
field equations that are modified by the presence of the spacetime defects. We would like to set 
up the modified field equations at a perturbative level. In the first part of the current section we 
will follow the lines of Ref. [9]. In the second part we will show that the results obtained can be 
reproduced with the help of the perturbative Feynman rules. These are given in App. [Bj where 



Eqs. (B.la) - (B.ld) can be directly derived from the action (2.1). The Feynman rule (B.le) 



follows from the deliberations of the previous section. 

A. Perturbative Ansatz for the photon field 



In order to obtain the modification of the photon field originating in the interaction with space- 
time defects via the scalar field <f> we have to solve the modified photon field equation resulting 
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from the action (2.1). In momentum space it can be written as an integral equation 

A (o) 



<?qg M {q)e^ qii (k - q) s A,^ a (k - q) . 



(4.1) 



\q^\>l/3l 
Ate{0,l,2,3} 



The index ',58" denotes that the system is at first considered in a finite rectangular region of side 
length St. The exact solution to the latter equation is out of reach. Therefore we make the following 
perturbative Ansatz for the full solution A u ^ in powers of the bare coupling constant \(°\ This 
is reasonable if we expect a modified photon dispersion law since current experimental bounds on 
Lorentz symmetry violation — and thus a modified dispersion relation for the photon — are very 
tight (see [T5]). 

A% = A^ v + \ {0) A^ U + {\ {0) ) 2 A ( § U + .... 

Herein A^ v is a solution of the free-field equation 

k 2 A^ u (k) = 0. 



(4.2) 



(4.3) 



By successively inserting the power expansion of Eq. (4.2) in Eq. (4.1), we obtain a perturbative 



expansion of the exact solution. Now the first step is to insert A^°\ We identify each perturbative 
order remembering that k 2 A^ u = 0. Using the definition 

A = — !— , (4.4) 
k 2 + le 

with an infinitesimal real parameter e to avoid the pole at k 2 = 0, we find the following first order 
perturbative solution, where we now take the limit M i— > oo 

A (0) 



#iW"(fc)= lim 



(2tt)' 



A(fc) / d 4 qg^q)e^q ll (k-q) g A^(q). 



(4.5) 



The first order photon field correction vanishes in the limit considered. This is clear from Eq. (3.21 ) 
since ]jsi^ ao {q) does not come together with a second background field, thus only producing a 
contribution for q = 0. 

The second order solution of the photon field equation reads as follows 

A (o) 



( A (0))2^(2), (fc) 



lim 



(2vr) 4 
( A (0))2 



A(fc) / A A qg^q)e^q,{k - q) Q (\<® A aM) 



lim ^-—^A(k) U\g Si {q)e^q ll {k-q) e A{k-q) 



(4.6) 



We contract the Levi-Civita tensors, perform the limit i— > oo, and use Eq. (3.21) to simplify the 
latter result which finally leads to 



(\^) 2 A^%k) = -C^A(k)B v Jk)A^(k) , C(°) = H(6(°)) 4 (A(°)) 2 ^, 



(4.7a) 



d A q 



1 2 



(b(°)y 



-Mi) 



1 k» 

(k-q) 2 + ie 7 



(4.7b) 
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with 



q) Q q a k? . 



(4.7c) 



In Eq. (4.7c), [a, 0, 7] denotes a totally antisymmetric permutation of the indices a, (3, and 7. The 
tensor K u also appears in the modified photon field that is obtained in the context of the effective 



background field model in |9|. The quantities 
hence is also a dimensionless parameter. 



i, (& (0) ) 4 p, and A(°> have no mass dimension and 



B. Establishing the perturbation from the Feynman rules of the modified theory 



Analogously, the perturbative series of the full solution of the modified photon field equation 
can be obtained in terms of Feynman diagrams. The corresponding Feynman rules follow from the 



action of Eq. (2.1 ) and are given in App. |B| 

We couple the photon field to a conserved fermion current j v (with k u j u (k) = 0) that is repre- 
sented by a plain line with an arrow, where the scalar field is denoted by a plain line without any 
arrow. The photon field is drawn as a single wiggly line and a double wiggly line stands for the full 
field. Ordinary vertices are represented by dots and the scattering at a defect ("defect vertex") is 
shown as a cross: 





V 


- + < 




k a 

1 Jv ' 


it 


jv ' 



v 



+ 

jv 




7 



A^°\k) 
+ ... 

it 



(4.8) 



Contributions proportional to odd powers of A^ ) (containing a defect vertex connected to only 
a single <f> field) vanish because of four-momentum conservation. For example, this is the case 
for the second diagram on the right-hand side of the diagrammatical equation above, cf. the 
discussion below Eq. (4.5) in the previous section. The second order perturbative solution A^ u {k) 
— corresponding to the third diagram on the right-hand side of Eq. ( 4.8 ) — leads to a nonvanishing 
correction of the photon field. 

A resummation of all one-particle reducible diagrams at one-loop order leads to a resummed 
photon field 

A"(k) e 



+ ... . 



(4.9) 



Multiplying the inverse standard photon propagator A _1 (fc) = k 2 with A^{k) nullifies the zeroth 
order contribution corresponding to the free photon field A^(°\k). Furthermore, it cancels a 
propagator in each further term and therefore an overall contribution can be factored out: 



+ 



(4.10) 
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This prefactor is the one-loop diagram in front of the round brackets on the right-hand side of 
Eq. ( |4.10[ ). It is given by the contribution -C^B^k) with B v ^{k) of Eq. ( |4.7b[ ). Hence, the 



modified field equation that results from the resummation of all one-loop photon field corrections 



is 



k 2 A u (k) = -C^B v Jk)A~<{k) . 



(4.11) 



A modified photon dispersion law following from this will be derived later. Note that one-particle 
irreducible higher-order corrections, for example 




are not covered by the resummed photon field A^. They are assumed to give a contribution to the 
modified dispersion relation of the photons as well. However since these are suppressed by at least 
one further factor (A^) 2 we will neglect them in our calculations. 

V. LEADING-ORDER PERTURBATION OF THE PHOTON FIELD 
A. Dimensional regularization 



We now want to compute the one-loop contribution B y [k) to the photon field that was set up 



in Eq. (4.7b). We begin by contracting the indices in Eq. (4.7c) 



K^A^{k) = SffiSfaik - q) K q a k^A^ 



2„2 



= k v {k ■ A^)q 2 - k v k e A^ a q eq(T - A^ v k~q 

+ A^ u k e k a q Qq(T - q v q e k s {k ■ A<®) + q v q Q A^ e k 2 . 

The second step is to perform the four-dimensional momentum integral over q: 



{2itfB%A^~i{k) 



(q 2 - 1 jb 2 + ie) 2 [(k - q) 2 + le] 7 
ga ^k s k a A^ u - k v k e A^ a + Tf a \k 2 A^ e - k e {k ■ 1^)1 } 



I. 
+ h 



k"(k-A^)-k 2 A^' 



where l Qa is a tensor one-loop and Iq a scalar one-loop integral 



I 



dV 



QgQa 



f-l/b 2 +ie) 2 [(k-q) 2 +ie] 



(5.1) 



(5.2) 



(5.3a) 



In = L 



d A q 



(q 2 -l/b 2 +ie) 2 [(k-q) 2 + ie] ' 



(5.3b) 



By power-counting we see that the integrals I ga and Iq are ultraviolet-divergent. Therefore they 
have to be regularized and we decide to use dimensional regularization. The basic principle is to 
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analytically continue the integrals to d spacetime dimensions, where d 7^ 4 is a real number. If we 
use the convention d = 4 — 2e, with four spacetime dimensions to be recovered in the limit e t— >• 0, 
the divergences become manifest as poles in e. Via 

/ d4q = (2 ^ )4 / wf " (27r) v " / (S 5 = (27r/i)4 " / drf(? ' (5 - 4) 

the renormalization scale of mass dimension 1, is introduced to conserve the dimension of the 
integral. 



B. Passarino Veltman decomposition 



Equation (5.3a) gives a tensor integral that can be reduced to scalar integrals I± and I 2 with 



the following Ansatz 

Iqo = rj ea Ii + k e k a I 2 . 
The integrals I\, I 2 follow from the contractions K\ = k e k° I ga i K 2 = ij^lga 

Y _Ki- k 2 K 2 y _ -dKi + k 2 K 2 
1_ {l-d)k 2 ' 2 ~ (1 - d)k* ' 



(5.5) 



(5.6) 



What remains is the reduction of the contractions K\ and K 2 to scalar master integrals via a 
Passarino- Veltman decomposition. This leads to the following result (see App. [C] for a detailed 
calculation) 



Ki = \{- A °{wf 



+ k z + 



1 



+ 



(6(0) ) 2 

2k 2 



-^[-k,—^-,o]-B (o, ' 



+ 



(6(o))2 (ft(o)) 



' (6(°)) 2 
C (-k,0, 



1 



0, 



(6(°)) 2 ' (6(°)) 2 
1 



(6(°)) 2 ' ' (6(°)) 



(5.7a) 



K 2 - /o - Bo { -k, -7^,0) + ttt^Cq ( —k, 0, 0, ' 



(6(°)) 2 ' 7 (6(°)) 2 u V ' ' (6 (0) ) 2 ' ' (b^) 2 , 
Here the contractions K±, K 2 are expressed solely in terms of master integrals given by 



B 0, 



A 
1 



1 



(6(°)) 2 
1 



(2vr^) 4 - d / d d q 



1 



q 2 - l/(6( ) ) 2 +ie' 

(2^/d^ (g2 _ 1/(b ; 0))2+ie)2 , 



(6(°)) 2 ' (6(°)) 2 

Bo (-^'°)=^ 4 "V d ^ 



1 



C °( fc ' ' ( 6 (0))2' ' ( 6 (0))I 



d 4 ^ 



/(6(°)) 2 + ie)p-g) 2 + ie] 
1 



(g 2 -l/(6(°)) 2 + ie) 2 p-g) 2 + ie] 



(5.7b) 



(5.8a) 
(5.8b) 
(5.8c) 
(5.8d) 
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Note that the Co-integral has neither infrared nor ultraviolet divergences so there is no need to reg- 
ularize it. The integrals can be computed by standard methods such as Feynman parametrization 
(see for example [16J). The results are as follows 



1 



1 



17T Z 



1 



1 

ITT 



2 B I 0, 



) 2 ' f 



(M°)) 2 

e ~ ln l(6(°))V 



In 
1 



(&(°))V 

+ 0(e) , 



+ 1 



+ 0{e) 



/o 
In 



dx In 
1 



fcV _ + l/(6(0))2) x + l/(&(0))2 _ i£ 



(5.9a) 
(5.9b) 

+ 0(e) 



(&(°))V 

fc 2 _ l/(fe(0))2 



+ 2 



A: 2 



In 



l_ (fo (°)) 2 A; 2 



le 



+ 0(e), 



1 

ITT 



2 C 



-jfe.O, 



r,0, 



(6(°)) 2 ' ' (6(o)) 5 



drr 



1 — x 



o /fc 2 x 2 - (/t 2 + l/(6(°)) 2 )x + l/(6(°)) s 



le 



k 2 



In 



1 - (6(°)) 2 A: 2 



le 



(5.9c) 



(5.9d) 



We have used 



-=~-1e + ln(4vr) , 



(5.10) 



with the Euler-Mascheroni constant 7 B ~ 0.577216, which is a reasonable redefinition of the regu- 
larization parameter. Terms of 0{e) have been discarded since they are not needed. An elaborate 
computation of the scalar integrals is presented in App. |Dj 

To summarize, we obtain the following photon field correction at second order in perturbation 
theory 

(2vfB^{k)A^\k) = ( Vea h + k e kj 2 ) 

x ik e k a A^ v -k v k e A^ a + rf a k 2 A^ e -k s (k-A^) } 



k»(k-A^)-k 2 A^' 



(k u k^ - r] u ^k 2 ) 



In - 2 Ji 



k 2 h 



A^(k) 



(5.11) 



The scalar integrals I±, I 2 result from the contractions K\, K2 via Eq. (5.6), where these are given 



by Eqs. (5.7a), (5.7b). The bare correction to the photon field is transverse and contains \je poles. 



In order to obtain a physically meaningful result Eq. (5.11) has to be renormalized. 



C. Renormalization procedure 

The second order solution of the photon field equation can now be written as follows 

(A(°))2A( 2 )"(fc) = -C (0) A(k)B^{k)A^(k) = -C {0) A(k)m^(k)A^\k), (5.12a) 

in^(fc) = i(jtT-ij^ 2 )n()t 2 ), (5.12b) 

H(k 2 ) = -i(/ - 2/i - k%) , (5.12c) 
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with the explicit result 



167r 2 n(/s 2 



1 

Ye 



2 



1 



(6( )^) 2 
2 



+ 



{b^kf 



(b(°)k) 2 (6(0) A;) 



In 



1 - (b^kf 



le 



(5.12d) 



Since the one-loop diagram computed resembles the vacuum polarization contribution of standard 
quantum electrodynamics (QED), we perform a renormalization of the coupling constant A^ '. We 
construct the modified photon propagator from the modified photon field by successively inserting 
one-loop corrections. This gives an infinite resummation of one-particle reducible diagrams at one- 
loop order leading to the full propagator. In this procedure we neglect higher order perturbative 
corrections 0((A(°)) 4 ) that are one-particle irreducible. Finally the photon propagator is coupled 
to a conserved current. This procedure will not be performed explicitly, but it is needed in order to 
drop all terms proportional to the four- momentum k^. In a diagrammatical notation our approach 
appears as follows 



+ 



(5.13) 



In terms of equations this corresponds to the vacuum expectation value of the time-ordered product 
of field operators (TA^(k)A u (k)), which is the resummed Feynman propagatoJ 4 ' 



{TA»{k)A v {k)) 



-it] 



k 2 



+ 



-a] 
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k 2 



+ 



-li] 



le 



ilV(fc) 



-irj u 



k 2 + ie 

act 



-17/" 



-17/ 



/3i/ 



li] 



k 2 + ie 'k 2 +ie apy 'k 2 + ie 

\W v {k) in^/^IF^A;) 



+ 



k 2 + ie k A + ie k e + ie 

The contraction of two transverse structures W i,u (k) results in 

U va (k)UJ(k) = -k 2 Yi vs . 



+ ... . 



(5.14) 



(5.15) 



So the resummation of all one-particle reducible diagrams at one-loop order corresponds to a 
geometric series and leads to the following result j^] 

(TA»(k)A v (k)) = {l - (A(°)) 2 II(fc 2 ) + [(X^) 2 U(k 2 )] 2 =f . . . j 

-irf 1 



fc 2 + iel + (A(°)) 2 n(fc 2 ) 



(5.16) 



The full propagator can be expressed via the bare propagator multiplied by a prefactor that contains 
the bare one-loop correction Yi(k 2 ) to the photon field. In order to perform the renormalization 



4 The coupling to the conserved current is suppressed. 

5 In the calculations performed within this section, (A^ ') 2 is understood to be extracted from the dimensionless 



constant C' ' defined in Eq. (4.7al. Consequently, it appears together with the one-loop correction il(fc 2 ) in order 
to keep track of all powers of A 
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we consider a physical process that contains a photon propagator, e.g. the scattering of <f> at a 
photon and its subsequent emission. We use the full propagator for setting up the amplitude of the 
process with the propagator momentum squared k 2 corresponding to the squared center of mass 
energy y/s, namely k 2 



s. 




(AW) 2 



-irj 



ft i' 



1 



s + ie 1 + (A(°)) 2 n(s) 



s. 



fll/ 



s + ie 



S, 



(5.17) 



Here S^ u contains the remainder of the amplitude, which is not important for current consider- 
ations. In the course of renormalization the bare coupling A^ ) is replaced by the renormalized 
coupling A, such that the renormalized amplitude is finite 

(A(°)) 2 A 2 



X 1 



(AW) 2 



A^ + 0(A 4 



(5.18) 



i + (A(°)) 2 n(s) ' K " ' i-A 2 n( s ) 

If the propagator momentum squared k 2 differs from the scale s we obtain an expression similar 
to the bare amplitude at order (A^ )) 2 



(AW) 5 



-177' 



fc 2 + iel + (A(°)) 2 n(A: 2 ) 



S, 



-IT]' 



flV 



k 2 + ie [l-A 2 n(s)] [l + X 2 U(k 2 )] 



Sau + 0(X 4 



-vf 1 



-s^ + o(x 4 



k 2 + iel + X 2 U mn (k 2 ) ^ v 
but with A^ ^ replaced by A and H(k 2 ) replaced by the renormalized one-loop correction 

n rcn (fc 2 ) = u(k 2 ) - u(s) . 



(5.19) 



(5.20) 



Equation (5.16) shows that the imaginary unit in front of the one-loop correction H(k 2 ) is put into 
the propagator. What in fact matters for the photon dispersion relation is the real quantity H(k 2 ). 



From the bare one-loop correction we obtain the renormalized correction according to Eq. (|5.20|): 

1 



16vr 2 n ren (fc 2 



2(6(°)) 2 



1 

+ 2 



+ 



1 



i 



(M°)fc) 2 (&(°)/fc) 4 

2 1 

(6(°)) 2 s + [(6(°)) 2 s] 2 



In 



1 - (b^kf 



ie 



In 



1 



(b^) 2 s-ie 



(5.21) 



where the renormalization scale fi obviously cancels. The physically important quantity H mn (k 2 ) 
is the difference between the bare self-energy correction evaluated at k 2 and the same quantity 
evaluated at an arbitrary scale s. In principle s can be chosen such that the last expression of 
Eq. (5.21) in rectangular brackets vanishes. This holds for s = l/(b^) 2 which then leads to 



16vr 2 n ren (A; 2 ) = 16vr 2 II ren (£; 2 ) 



1 



1 



3=l/(6(°)) 2 

I ^ [l_(6(0) fc )2] 



In 



1 - (b^kf 



ie 



(5.22) 



2 V" (b(°)k) 2 J 2(&(°)£;) 4 
The choice s = 1/(M°)) 2 is not unreasonable since l/(&(°)) 2 is the only parameter of the action 



(2.1) that has the same mass dimension as s. 
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VI. LEADING ORDER PERTURBATION OF THE SCALAR FIELD 



Having obtained the modified photon field we now consider the modification of the scalar field 



resulting from the action (2.1). There are two corrections to the bare scalar field. First of all 



there is a quantum correction involving the photon (see Fig. 2(a) ) and, secondly, the scattering of 



the scalar field with the defects has to be taken into account (see Fig. 2(b) ). We will now compute 
the self-energy correction given by the diagram in Fig. |2(a)| 



-i£(fc 2 ,m^ 




(iA) 2 (27r/i) 4 - d J d d qe a ^q a (-k - q) p e^(-q) y (k + q)g 

x -hy -h/ gCT 

q 2 + m 2 + ie (k + q) 2 + m 2 + ie 

2X 2 [k%(m) -T 4 (m,k)} , 



(6.1a) 



with 



I 3 (k 2 ,m 2 ' 



(2irrf- d J d d q 



(q 2 — m 2 + ie)[(/c + q) 2 — m 2 + ie] ' 



(6.1b) 



h(k 2 ,m 2 ) = (2tt//) 



4-d 



d d q 



(k ■ qf 



(q 2 — m 2 + ie)[(A; + q) 2 — m 2 + ie] 



(6.1c) 



The ultraviolet divergences are again regularized by dimensional regularization with renormaliza- 
tion scale [i. Furthermore a photon mass m has been introduced to regularize possible infrared 
divergences. The results of the integrals ^3, I4 in the limit m 1— > are given by 



I 3 {k 2 ,0) = 0, h(k 2 ,0)=m'- 



k 4 ri 



+ 



In 



ie 



(6.2) 



Fortunately the integrals are infrared-finite. Hence we consider T,(k 2 ,0) = S(/c 2 ) from now on. A 
resummation of all one-loop photon self-energy corrections leads to a modification of the scalar 




(a) One-loop photon 
self-energy contribution to 
the scalar field. 



(b) Correction to the 
scalar field originating 
from the scattering of (j> 
at the defects. 



FIG. 2. Possible corrections to the bare scalar field (j)(k). 
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field propagator 



\b 2 



k 2 - l/(6(°)) 2 - (6(°)) 2 S(P, 0) + ie 



+ •- 



— o 
oo 



-•+ 



(6.3) 



Contrary to Sec. [V] the renormalization of the correction to the scalar field will be performed with 
the help of counterterms, which is a more convenient procedure here. We will follow the lines of 
|17j . Both the photon mass and the field are renormalized according to 

/ — Z 2 1 

= V Z 2 (pren , h/b = JU^y ~ ' $Z 2 = Z 2 ~ I , (6-4) 

where <j) is the bare scalar field, (p ren the renormalized field, Z 2 the field renormalization constant, 
5z 2 the field renormalization counterterm, and 5i/i the mass counterterm. Furthermore, we use 
the renormalization conditions 

d 



)|fe2=l/ 6 2 = 0. 



d/fc 2 



S(r) =0, (6.5) 

fe 2 =l/fe 2 

From these we obtain the counterterm which can be used for both the mass and the field renor- 
malization 

'k 4 



with 



h/b 



h/b 2 



$z 2 



7T 



7T 



2A 2 



+ A 2 



2A 2 



A 2 

V + T 



2 In 



13 - 8 In 



1 



6 2 /i 2 
1 



ie 



ie 



(6.6a) 

(6.6b) 
(6.6c) 



The renormalized photon one-loop self-energy contribution to the scalar field is then 

'k 4 



-iE(r) + i 
ivr 2 A 2 . j 



46 4 



6 4 / ft 4 



l + 21n 



1 



+ ie 



(6.7) 



This correction leads to a renormalization of the 0- field mass. The correction is of order A 2 and 
depends on both k 2 and the renormalized mass 1/6 2 



1 



1 



^ ^ + b^ Ten (k 2 ,b 2 ). 



(6.8) 



(6(0) )2 6 2 

As in the case of the photon field we now define the resummation of all one-particle reducible 
one-loop corrections to the scalar field 



+ 



O CM>- 



+ •- 



+ counterterms + . . . . 



+ 



(6.9) 
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This resummation fulfills a modified field equation with renormalized mass 1/6 of the eft- field 



1 

u 2 



l 

V 2 



4>(k) 



O 



d>(k) 



(6.10) 



Inserting this mass correction into the photon field equation leads to 0(A 4 )-corrections. Note that 



the imaginary part resulting from both Eq. (6.1a) and Eq. (6.7) is given by 



Im(E) 



7T" 



2 lA 



-X z k 



(6.11) 



This result does not depend on the renormalization program since by the optical theorem it is 
linked to the total cross section for the decay of an excitation of the (/>-field (a scalar particle) into 
two photons. 

Finally, we take into account the second contribution in Fig. 2(b)[ namely the scattering of <f> 
at a defect. In so doing we define the full solution 4> (k) of the corresponding field equation as the 
resummation of all one-loop corrections with defect vertex insertions 



• +« X +— X- 



O O 



+ 



+•- 



+ 



+ counterterms + . . . . 
This leads to the following field equation for the scalar field 



b 2 V b 2 



(6.12) 



(6.13a) 



or 



4> s (k) = Q(p s (k) . 



(6.13b) 



b 2 v b\ 

Physically, the interaction of the (/>-field with the defects results in a shift of the mass of (j> that 
corresponds to the density of defects as long as momentum transfer to the defects is neglected. 

' (6-14) 



1 1 L 2 



b 2 

Q 



VII. MODIFIED PHOTON DISPERSION RELATION 



A. Final result for the modified theory 

To summarize, the second order perturbative solution of the photon field equation is 
X 2 A^ u (k) = -CA(k)B^(k)A^(k) = -CA(k)iU^(k)A^{k) , 

mr{k) = i(k u k^ -r]^k 2 )u rcn (k 2 ), 



16vr 2 II ren (fc 2 



{b e kf 



[i - (MO 

2{b e kY 



21 2 



In 



1 - {b e kf 



le 



(7.1a) 
(7.1b) 

(7.1c) 
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where C = Eb 4 g X 2 g and e = 0+. In Sec.M we integrated out the scalar field at one-loop level. Using 
the results from Sec. VI this leads to an effective vertex whose Feynman rule is given by Eq. (B.lf) 
in App. [Bj 

In order to obtain the modified photon field equation, we define the one-loop resummed photon 
field A v e (k) by employing the full scalar field solution 4> g (k) from the previous section. 



Am) 



~t 



+ 



(7.2) 



k 



k 



k 



k 



k 



k 



The one-loop contribution to the photon field is transverse, and after renormalizing the coupling 
constant using Eq. (5.19) it is finite. Furthermore, it becomes imaginary for k 2 = 1/b 2 . This 



indicates that an electromagnetic wave is damped when k 2 approaches the mass of the scalar field 
4> leading to a resonance behavior. 



Pursuing the discussions at the end of Sec. IV B and below Eq. (5.20) we now consider the 
modified field equation 



^^^^ 



(7.3a) 



k 2 A u g {k) = -C[- W (k)]A20e) = -C (k'ky - 5%k 2 )n rcn (k 2 )Al(k) 



Ck 2 % en {k 2 )A v Ak). 



(7.3b) 



Two distinct photon dispersion relations follow from Eq. (7.3b). First of all we obtain the standard 
dispersion law k 2 = 0. Secondly, considering k 2 ^ Eq. (7.3b) results in 



k 



A v e {k) = -CA(k)[- iB" (k)]A2(k) = C-^—fl ren (k 2 )A»(k) . 



(7.4) 



k 2 + ie 

Nontrivial solutions for the photon field will exist if the following transcendental equation holds 

i-cn rc # 2 ) = o. (7.5) 

Unfortunately it cannot be solved analytically. Since we expect the solution to be a minuscule 
correction to the standard photon dispersion relation k 2 = 0, we assume (b g k) 2 < 1. This leads to 
an approximate equation that can be investigated analytically: 



1 + 



C 

64vr 2 



2b 2 

1-f (*o 2 -|k| 2 ) 



0. 



(7.6) 



Its solution is 



k| 2 + rriy , 



~C~ 



b 2 

Q 



(7.7) 



with the three-momentum k = (^1,^2,^3). Thus the photon acquires a mass m Y . Note that this 
mass depends inversely on the parameters b e and C, which are supposed to be small. 
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B. Physical meaning of the modified photon dispersion relation 



In this section we wish to discuss what the result obtained in the previous section means 
physically. For our argument we choose exemplary values for the parameters of the theory. First 
of all, a reasonable value for the upper boost limit A' has to be chosen. This is a dimensionless 
number, thus e.g. the ratio of two characteristic masses. As the first we choose the Planck mass 
Mpi = ^hc/G w 1.22 • 10 19 GeV/c 2 and as the second the proton mass M P « 0.94GeV/c 2 that is 
a characteristic scale for standard model physics. Since the Lorentz group acts on four-dimensional 
Minkowski space it can be shown that H grows with the fourth power of the dimensionless number 
Mpi/Mp that has been previously mentioned. 

Furthermore, since the scalar field is the mediator between the photon and the spacetime defects, 
which are Planck-scale effects, the scalar field mass l/b e is assumed to be large, perhaps some 
fraction of the Planck mass. Then b e would be of the order of the Planck length. Assuming the 
spacetime defects to have an average separation of 10 10 x Lp\ (see the discussion below Eq. ( 3. 13| ) ) 
and setting A ~ 10~ 19 (the coupling constant A is taken to be a very small number) we obtain: 



C 



IgbtX 2 ~ 



1 



100 \M Pl /M P J VV(10 10 x L P1 ) 4 




(7.8) 



Hence, the upper choice of the values is consistent with our procedure of performing a perturbative 
expansion with respect to A, which in principle corresponds to an expansion in C. 



Now it becomes evident that the behavior of Eq. (7.7) as a photon mass is unusual for a 



perturbative approach. The deviation from the standard photon dispersion law is expected to be 
a correction that vanishes for b e \— > or C i— > 0. Physically, for 6„ i — >• the scalar mediator between 
the photon and the defects decouples, whereas for C t- > either the density g of defects or the 
coupling A between the photon and the scalar vanishes (or again b e t— >• 0). For both scenarios the 
modified photon mass is expected to vanish, since it is a perturbative result. However, as becomes 



evident from Eq. (7.7), the contrary happens. 



Let us see how this behavior is connected to the limit A; 2 i — > of the quantity n ren of Eq. (7.4), 
which is the mathematical basis of the modified photon dispersion law. Two cases have to be 
considered here: {b g k) 2 -C 1 and (b g k) 2 S> 1. 



lim 16vr 2 k 2 U reTL (k 2 



(fe e fc) 2 <i 

fixed 



o, 



(7.9a) 
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+ lim —z I 1 

A: 2 1-^0 b 2 g 



2b\k 2 



ln(l - ie) 



(7.9b) 



Compare these results to the renormalized self-energy correction of the photon (vacuum polariza- 
tion) n ren (/c 2 )| < ^ ED in ordinary QED, for which 



lim fc 2 n ren (/c 2 ) 
k 2 i-4-O 



QED 



0. 



(7.10) 



Equation (7.10) means that gauge invariance is maintained for quantum corrections in QED and 
the dispersion law of the photon remains as k 2 = 0. An analogous argument holds for the modified 
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theory defined by the action (2.4). As long as k 2 -C 1/& 2 ,, for which the scalar mass is large and 
the modification of the photon dispersion law is a small deviation from k 2 = 0, gauge invariance 
is maintained. Then the photon mass given by Eq. (7.7) cannot appear in this case. The tensor 
structure (k v k^ — r] u ^k 2 ) of B l/ ' Y (k) in Eq. (7.1) is crucial for this result, namely the conservation 
of gauge invariance. As long as this structure is conserved by the interaction with the defects, the 
photon dispersion relation stays k 2 = (for b 2 k 2 <C 1). 

However from Eq. (7.9b) it follows that gauge invariance is violated for large b e (and small 
scalar mass l/b g ). The technical reason is that e = + cannot be discarded in this case, since the 
argument of the logarithm in Eq. (7.1) may be negative. The complex logarithm has a branch 
cut on the negative real axis and, therefore, a small imaginary part e = + has to be added to its 
argument. Thus the limit k 2 H > does not exist here and gauge invariance is violated resulting 
in a photon mass. Because of the infinitesimal imaginary part in the logarithm, n ren (fc 2 ) also has 
an imaginary part. Physically this corresponds to the damping of electromagnetic waves when the 
modified photon momentum square approaches the mass square of the scalar field. 



C. General remarks about the previous results 



A graphical investigation of n ren (fc 2 ) shows that Eq. (7.5) does indeed not have a solution. This 



means that the perturbative solution (17. 71) only appears as a solution to the expanded equation 



(7.6), but it is not a solution of Eq. (7.5), which is exact. Why did we not do a graphical analysis 
at first? The answer is that in this section we would like to argue on a very fundamental basis 
why the photon dispersion relation remains the conventional one within the simple spacetime foam 
model proposed. For this reason we will investigate a generic function h coming from a quantum 
correction that involves the interaction with the defects. 

However first of all we recall the following three general modifications of the photon dispersion 
law that originate from a violation of Lorentz or gauge invariance. They are denoted as cases (1), 
(2), and (3). 

1) Scaleless Lorentz-violating modification (cf. modified Maxwell theory |18H21j ): 

In the first case Lorentz invariance is violated by preferred spacetime directions. These are 
given by four- vectors a%, a 2 etc. with dimensionless components «i, k 2 etc. The modified 
photon dispersion relation is obtained from a zero of a function / whose arguments are 
scalars. Therefore, they can only consist of scalar products with the momentum four-vector 
and scalar products among themselves: f(k 2 , a± ■ k, a 2 • k, a% • a 2 , • . . ) = 0. The solution of 
such an equation is 

ko = g(h,k 2 , k 3 , m, k 2 , . . . ) ~ |k| + 9i(h,k 2 , k 3 )Ki + g 2 (ki, k 2 , k 3 )n 2 + ... , (7.11) 

with functions g±, g 2 etc. that have mass dimension 1 and only depend on the three- 
momentum components. Higher-order terms in the scaleless deformation parameters have 
been neglected. 

2) Scale-dependent Lorentz-violating modification (cf. Maxwell-Chern-Simons theory |19H22j ): 

Preferred spacetime directions a\, a 2 , etc. appear in this case as well. They come together 
with physical scales mi, m 2 etc., hence have dimensionful components. We will denote them 
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with Latin letters in the following. The modified dispersion relation is a zero of a function 
/ containing only scalar products as arguments: f{k 2 ,a\ ■ k, a 2 • k, a± ■ a 2 , ■ ■ ■ ) = 0. We 
perform an expansion in the scales mi, m 2 of the deformation by considering them to be 
much smaller than the remaining physical scales. This leads to the following approximate 
solution of the latter equation: 

ko = g(h,k 2 , k 3 , mi, m2, ■ ■ ■) ~ |k| + <7i(fci, k 2 ,h)mi + g 2 (h, k 2 , k^)m 2 + ... , (7.12) 

where the functions g\, g 2 etc. are dimensionless and depend on the three-momentum com- 
ponents. Thus the existence of preferred directions leads to a different behavior for the 
four- momentum components in the dispersion relation since the scalar products project out 
certain momentum components. 

3) Emergence of a photon mass: 

The emergence of a photon mass m y by a violation of gauge invariance would result in the 
following dispersion relation 

2 

1 in 

k = ^/k 2 + m 2 = |k| + — ^ + ... . (7.13) 



Note that in Eq. (7.13) no term that is proportional to m y appears — contrary to Eq. (7.12) 



where linear terms in the dimensionful parameters mi, m 2 etc. can be found. 

A modification of the photon dispersion relation in the context of spacetime foam is expected to 
emerge according to one of the previous cases. For example, within the spacetime foam model 
considered in [9] the photon dispersion law is modified according to (1). The reason is that the 
method of distributing spacetime defects in the latter reference is not Lorentz- invariant, leading to 
the spatial dimensions behaving differently to the time dimension. This corresponds to a preferred 



direction in spacetime and results in a modified photon dispersion law of the form of Eq. (7.11), 
since in the photon model there is no quantity having a mass dimension. 

On the contrary, cases (1) and (2) cannot play a role within the spacetime foam considered in the 



current article. In the action (2.4) of the effective theory no preferred spacetime directions appear. 
Furthermore, no such directions emerge, since the defects are distributed in a Lorentz-invariant 
way. As a result of this time and spatial dimensions behave the same. Thus only the dimensionless 
quantity b g k = b Q (kQ — k 2 ) 1 / 2 appears in the equation leading to the modified photon dispersion 
law. So it must be of the form h(b g k) = 0, where h is contains the expression of the quantum 
correction involving a scattering with the defects. In our case h corresponds to the left-hand side 
of Eq. (7.5). Assuming b 2 k 2 <C 1, which is reasonable within perturbation theory, h(b g k) can be 
expanded with respect to its argument b e k to quadratic order: 

MO) + b e kh'(o) + ifiJfcVco) ~ =► k 2 ~ k 2 - , (7.14) 

where we have used that h'(0) = oj^] This would lead to a nonvanishing photon mass indicating a 
violation of gauge invariance. Hence, the only modification may be according to (3). However the 



6 Avoiding branch cuts by infinitesimal imaginary parts at the appropriate places, perturbative corrections are 
analytic functions in k 2 , and so do not depend on k. 
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photon mass in Eq. (7.14) is proportional to the scalar mass l/b g , which is not a small perturbation 
when l/bg is large. Considering the left-hand side of Eq. (7.5) as the function h we, in fact, obtain 
the result of Eq. (7.7). Hence even without a graphical analysis, which is not possible for a general 
function h, the argument given is still valid. We conjecture that in such a case, where the photon 
interacts with the defects via a quantum correction with a scalar particle, its dispersion relation 
remains conventional. 

What is the physical reason ruling out a modification (3)? Because of £; 2 n ren (&; 2 ) = for A; 2 i — > 



(see Eq. (7.9a) gauge invariance is indeed conserved, although the action ( |2.1[ ) seemed to violate 
gauge invariance (with respect to the gauge transformation <fi \— )■ —(f)) explicitly. However, the 
sign change of the penultimate term of (2.1) can be absorbed into the defect charges e. For an 
infinite number of defects this does not change anything since there are equally many defects with 
e = 1 and e = — 1. Also for the last term in the action the sign can be absorbed into the coupling 
constant A. Since the correction computed is proportional to A 2 this sign change has no physical 
implications. 



VIII. MOMENTUM TRANSFER FROM PARTICLES TO DEFECTS 



The final result of Sec. VII was obtained after making a series of assumptions in the course of 



the calculation. Let us recap these assumptions: 

• Ass. (1) : The interaction of the photon with the defects is mediated by a scalar field according 



to the effective theory of Eq. (2.4). This forms the basis of the simple model proposed in 



this article and is, therefore, the most important assumption. 

Ass. (2a) and (2b) + (2c) + (2d): Isotropic and homogenous ("random") defect distribution 



(Sec. |III A) + dense defect distribution (Sec. IIIB) + infinite spacetime volume (Sec. IIIB). 



These have been introduced not only in order to keep the calculation feasible, but also for 
physical reasons. In principle one or several of them can be dropped leading to a more 
difficult calculation. 

• Ass. (3): The modified photon momentum squared is supposed to be much smaller than the 
mass of the scalar field squared, i.e. k 2 <C 1/b 2 . This assumption has been made only for 
physical reasons. The mass of the postulated scalar field l/b e is expected to be large and 
the deviation of the photon momentum squared from the standard result k 2 = should be 
a small perturbation. 

So far the momentum transfer from the photon to the defects and vice versa has been neglected. 



The reason for this is Ass. (2) and Eq. (3.21) that followed from it. 

Now we would like to investigate how the results obtained change when momentum is transferred 
from and to the defects. This means that we drop Ass. (2a), (2b), (2c) or (2d). The technical 



problem is that the free photon field in Eq. (4.6) then depends on the integration momentum p, 



which renders the evaluation of the corresponding integral impossible. However we will stick with 
Ass. (1) and treat the propagation of photons through a spacetime foam via the effective theory 



defined by Eq. (2.1). Thus the external photon momentum is assumed to be much smaller than 
the Planck scale. 
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On the one hand, the (/>- field with its renormalized mass 1/6 directly interacts with the spacetime 
defects. We assume that (ft will only probe the defects if 1/6 is of the order of the Planck scale. For 
1/6 much smaller than the Planck scale a momentum transfer to the defects will be low at most, 
i.e. suppressed by that scale. As a result, the main contribution of the p-integral will come from 
the region p ~ k. The difference between p and k must be suppressed by a small dimensionless 
number that can be written as a ratio of two mass scales. The scale where the influence of any 
spacetime foam may become especially important is the Planck mass Mp\. Since the (ft- field is 
assumed to interact with the defects, the only other scale is the mass 1/6. 

On the other hand, the interaction of the photon with the spacetime defects is mediated via 
a quantum correction involving the virtual (ft- field with the new mass l/b e originating from the 
interaction of (ft with the defects. As a result, the photon acquires a size from this quantum 
correction that is inversely proportional to l/b g . Hence we assume the suppression of a momentum 
transfer from the photon to the defects to be similar as for the scalar field, but with the mass 1/6 
replaced by l/b g . This behavior will be summarized in the following paragraph. 

A. Assumption (4): Momentum transfer suppressed by the Planck scale 



In order to be able to compute the integral over p in Eq. (4.6 ), we introduce a fourth assumption. 
We assume that the ratios 

(fc - p)2 ~6^p? {k - p) - X ~bki k - X > (8 ' 1} 
approximately hold for the momentum transfer k— p of a low-energy 0-field with initial momentum 
k and final momentum p to a spacetime defect. Analogous ratios hold for the momentum transfer 
of a photon to a defect with the difference that 1/6 must be replaced by l/b Q 

(t-rt—j^*-*. (8.2) 



In Eqs. (8.1), (8.2) x is an arbitrary but fixed four- vector in configuration space. 



B. Defect distribution with large separation between individual defects 



In Eq. (3.13) we assumed the distribution of spacetime defects to be dense (see the left panel of 
Fig.[3j). In this case the defect distribution can be approximated by an effective background field 
and the sum over all defects results in a spatial integral. Discarding Ass. (2c) means that the sum 
mentioned can no longer be approximated by such an integral and we cannot define a background 



field, any more. As a result, the 5-function in Eq. (3.13), which helps to get rid of the second 



integral in Eq. (4.6), has to be replaced by a sum again. However, due to Ass. (4) the momentum 
transfer is suppressed by the Planck scale and the integration momentum p does not appear 

J d 4 pg5(k-p)f(p)^^exp(-i^-k-x^f(k). (8.3) 

2 = 1 



The mass of the 0-field given by Eq. (6.14) then changes as follows 



1 1 b 2 ^ ( 1 , \ 1 , . 



b 2 b 2 y Z_^~ *- ^ bMpl ~*J ~ b ( k>x 
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FIG. 3. Each panel in the current figure shows a distribution of spacetime defects. The scales of the panels 
are chosen in order to make clear that the distances between the defects in the left square correspond to 
multiples of the Planck scale, whereas in the right square they are macroscopic. Hence the right panel 
illustrates a foam in which individual defects are separated by distances that lie many orders of magnitude 
above the typical distances in the left panel. 



It now explicitly depends on the ^-momentum k and the defect positions Xi (see the right panel of 



Fig. In the loop integral of Eq. (4.7b) the shifted mass has to be inserted. Then the integral 



depends on the defect positions Xj. This can be stated as follows 



B\(k)^B^(k,x l ) = Y t J^^P 



1 



b(q,Xi)M Pl 
1 



" Iv ' X- 



1 2 



b{q,Xi)< 



H(q,Xi) 



1 K* 

(k-q) 2 +ie 7 



with 



b(q,Xi) 



jH(q,Xi 



q 2 - l/b 2 (q,Xi) + ie 



and K% of Eq. (4.7c). 



.5a) 



3.5b) 



So the result of the loop integral B u " / (k) will also depend on the positions x\ of the defects. 
Because of this its tensor structure may involve terms such ^k 1 etc. and they will spoil 

the form of the standard tensor structure (k v k^ — rj u ^k 2 ). This leads to a dispersion relation of the 
photon that differs from k 2 = 0. 



C. Defect distribution in a finite spacetime volume 

Now we stick to Ass. (2c) but drop Ass. (2d) and consider defects in a spacetime with finite 
volume V. For computational reasons the shape of the volume shall be a four-dimensional cubd^] 



7 with side length £$, and centered at the origin 
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containing J\f defects. Because of the finiteness of V the positions Xi are bounded and the argument 
of the exponential function can be assumed to be small. So we expand the complex exponential 
function in Eq. (8.3) with respect to the (small) momentum transfer. In the remainder f(p) of the 



integral over p in Eq. (4.6) we effectively replace all p by k and drop the integral over p. 



This procedure leads to a new replacement rule for the 5-function 
d 4 p g5(k - p)f(p) ^±[M + jrS n 



V 



f(k), 



n=l 



.6a) 



3.6b) 



i=l 



Since we assume the treatment of the spacetime defect distribution to be a good approximation 



the sum in Eq. (8.6b) can be replaced by integrals. These are restricted to the finite spacetime 
volume V. 



J d 4 Pe 5(k-p)f(p) ^ 6 (l + 



In 



. x r 



.7a) 



3.7b) 



n\MV j v 
It is evident that this expansion only makes sense when the integrals run over a finite space, 
otherwise the integrals would be divergent as the integrands are not suppressed for x \— > oo. 

If we also keep Ass. (2a) and (2b) for the moment we end up with a "random distribution" of 
spacetime defects. This means that no preferred direction is defined by the distribution. In this 
case all integrals that are odd in x vanish, I-ik-x = f° r k £ {1, 2, . . . }, since there is no counterpart 
that could make up the index structure of the result. 

The tensor structure of the even integrals can be generated by the metric tensor only. That is 
why the even integrals do not vanish but depend on combinations of metric tensors. For example, 
in four spacetime dimensions the integral with two indices is given by 



d xx 2 



v -@ 2 V. 
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Neglecting constant prefactors, the behavior of integrals that are even in x is as follows 
So further terms in the expansion are suppressed by powers of ratios of the Planck scale and the 



2k 



ke {1,2,...}. 



(8.9) 



square root of the spacetime volume. Finally Eq. (6.14) has to be replaced by 

k 



1 1 u2 b 2 ^^ ( W V 

^^ 2 + v^>* ^PmpJ 



1 



3.10) 



where C<ik are mere numbers. Since the new mass of 4> neither involves any preferred directions 



nor any defect positions x% the resulting B V1 of Eq. (4.7b) will still have the gauge-invariant tensor 
structure {k h 'k 1 — ff^ 'k 2 ) . Although the shape of a finite spacetime volume is not Lorentz-invariant, 
the photon is not affected by the finiteness of the spacetime. The reason is that we do not set 
any boundary conditions on the photon field in the framework of the simple spacetime foam model 
considered. Thus the dispersion relation of the photon stays k 2 = 0. 
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D. Anisotropic or inhomogeneous distribution 

If we additionally drop Ass. (2a), the spacetime defect distribution may be anisotropic and, 
therefore, define preferred directions in spacetime (see the first two panels of Fig. [4]). Let us 
assume that there is one such direction: (£ M ) = £ , £ 2 , £ 3 ). Then the result of all spatial 
integrals will involve £. For example the result for the integral with two indices is then made up 
of the metric tensor and the tensor product of preferred directions 



/ 



d 4 xx^ = -^(cf (OT + Cf (OCT) * US.! la) 




(8.11b) 



(a) one preferred spacetime direction (b) two preferred spacetime 

directions 
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(c) inhomogeneous distribution 



FIG. 4. Here, the first panel illustrates an anisotropic distribution defining a single preferred direction £, 
whereas the second panel contains a regular lattice of defects defining two preferred directions £i and ^2- 
The third panel depicts a section of an inhomogeneous defect distribution, where in the center of the region 
shown the density is higher than near the margin. 
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Using Ass. (4), the second integral in the expansion can be written as follows 



h ~ ~W (^ 0) (0(*-p) a + cf\()[{k- P ) ■ c?) 



1 w 



PI 



U2) 



M P1 

It can be shown that despite the existence of a preferred spacetime direction the integrals i2fc+i still 
vanish because they are odd with respect to the integration variable. Since the tensor structure 
of each even integral can involve at the most 2k preferred directions £, we obtain the following 
general behavior for the new mass of the </>-field 

1,2 



^E(^)E42(0(*-<f'(^)^, (s.i3, 

where C^iC) are functions with respect to the preferred direction The loop integral now depends 
on £ and its tensor structure may involve terms such as C^C 7 ) Ck" 1 etc., which destroy the standard 
structure (k v k^ — ij^k 2 ). This case may also produce a deviation from k 2 = in the dispersion 
relation of the photon. For example, the appearance of a preferred spacetime direction is the reason 
for the modification computed in [9]. 

Although the calculation presented here has been performed for a finite volume there is nothing 
to suggest that the physical effect originates from this fact. We have just seen that a finite volume 
does not have any effect on the dispersion relation of the photon so long as the defect distribution 
is isotropic and homogeneous. 

Keeping Ass. (2a), but dropping (2b) would result in an inhomogeneous defect distribution (see 
the third panel of Fig. Hi). This means that the defect density cannot be considered as constant 



and, therefore, it cannot be pulled in front of the integral in Eq. (3.13). There is no physical input 



for a varying density function g = g(x), since we are not aware of a mechanism resulting in an 



inhomogeneous distribution of defects. The resulting integral in Eq. (3.13) will certainly be much 
more complicated if g is a function of x. If there are regions with a high (constant) density g^ 
and regions with a low (constant) density gi such that g^ ^> gi, the main contributions from the 
integral will come from the regions with g{x) = g^. Effectively, this means that the integration 
volume is reduced. How the photon dispersion relation is affected by such changes, will not be 
further investigated in this paper. 

IX. PT-SYMMETRIC EXTENSION AND THE PERCOLATION OF DEFECTS 

In this section we would like to deliver a brief discussion on an interesting issue. Replacing the 
real coupling constant A^ ^ in the action dTTl) by an imaginary one we find that the equation 



A: 2 [l + i 2 Cn ren (A: 2 )] =0 (9.1) 
has a second physically acceptable solution different from the standard dispersion relation: 

k 2 = a{ 1 )^, (9.2) 
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TABLE III. Transformation properties of the vector potential A^, field strength tensor F^", dual field 
strength tensor F^ v , four-momentum fc^, and the imaginary unit with respect to P, C, T, and their combi- 
nations PT and CPT. 




FIG. 5. Several data points obtained for the coefficient 0(7) in Eq. (9.2) are compared with the expected 
behavior for the order parameter: a(j) oc (1 — J c /j) ■ The exponent used is f3 = 0.64 which is also obtained 
in the context of a four-dimensional lattice percolation. From the plot it becomes evident that j c ~ 2. The 
latter point depends on the renormalization scale we choose. 



with a function 0(7) and 7 oc A 2 £. Equation (9.2) means that the photon becomes massive. 

This solution will exist if 7 is larger than a critical value j c . Furthermore a seems to lie in 
the interval [0,1]. The replacement — > iA^ makes the interacting part of the Lagrangian 
non-Hermitian. However this may not be a problem so long as the Lagrangian is symmetric under 
a combined parity transformation P and time reversal transformation T: (PT)C = C. Non- 
Hermitian but PT-symmetric quantum mechanics has been thoroughly studied over the past few 
years [23]. It was found that it can serve as a description of real physical systems, see e.g. [23H26] . 
PT-symmetric quantum field theories have not been investigated that profoundly. Nevertheless 
there are indications that such extensions of ordinary Hermitian theories are physically mean- 
ingful |27j . According to Table III the last term of the action (2.1) is indeed PT-symmetric for 
imaginary coupling constant A^. 

Observing the relation between 7 and g we speculate that this solution could be related to 
the percolation of defects. Percolation theory studies the formation and properties of clusters of 
objects randomly distributed over a lattice or a continuous space. This theory exhibits a phase 
transition. As the density of the objects increases so does the mean size of the clusters, and it 
becomes infinitely large at some critical density g c . For densities smaller than g c only clusters 
of finite size exist, whilst for densities larger than g c clusters of infinite size also appear. At the 
critical point the system is scale invariant. 

This assumption is supported by the observation that the behavior of 0(7) is well described by 
the critical exponent /3 of the percolation phase transition in 4 dimensions (cf. Fig. [5]). This issue 
deserves further study. 
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Here we see that the general argument about the occurrence of a photon mass given at the end 
of Sec. (VIIC) does not hold. This is because the argument relied on a photon mass continuously 
depending on the model parameters, which is not the case here (as can be seen in Fig. foh. 



X. SUMMARY AND CONCLUSIONS 

To summarize, we have investigated photon propagation through a spacetime foam made up of 



a Lorentz-invariant distribution of time-dependent, pointlike defects. From the action of Eq. (2.4), 
photons do not interact directly with the defects but via a scalar mediator field. 

If the spacetime volume is assumed to be infinite and the distribution of defects is dense, 
isotropic, and homogeneous, the only effect on the propagation of the scalar mediator is that 
of enlarging its mass. The defects then effectively act as a modified background analogous to 
the pictorial representation of a particle getting its mass via the interaction with the Higgs field 
vacuum. So for the scalar field the defects have been "integrated out" and the photon interacts 
with a scalar that has a mass larger than when compared to the case without defects. 

The surprising result is that the photon dispersion law remains standard: k 2 = k 2 . — |k| 2 = 0. 
Hence the photon does not "feel" the random background field that mimics the defects as an 
effective theory — at least at leading order in the interaction between the photon and the scalar 
field. The reason for this is that on the one hand no preferred spacetime directions will appear if 
the defects are distributed in a Lorentz-invariant manner. On the other hand the photon does not 
become massive as gauge invariance is not violated. In order to rule out the appearance of a gauge 
anomaly the perturbative correction to the photon dispersion law must indeed be computed, which 
leads to the nontrivial result obtained. 

The outcome of the calculation sheds some new light on the very restrictive bounds on Lorentz- 
violating parameters of the Standard Model Extension [15J. If a nontrivial spacetime structure 
is assumed to be the underlying cause for a possible Lorentz violation, low-energy experiments 
are unlikely to detect Lorentz violation if at energies much smaller than the Planck energy this 
structure can be described by the effective theory considered. 

However one also has to keep in mind that the photon dispersion relation continues to be the 
conventional one because of several idealized assumptions taken to ensure a feasible computation. 
Using physical principles we have tried to predict the effect on the photon dispersion relation when 
one of the assumptions is dropped. Discarding the assumption of . . . 



. . . the action of Eq. (2.4): The model is no longer an appropriate description as the defect 
structure itself is expected to become important. Furthermore, the theory could be modified 
such that the scalar field is discarded altogether and the photon be made to interact directly 
with the defects. Then it would probe the defects with its wavelength. If we still assume the 
momentum transfer (A/c) 2 to the defects to be suppressed by one power of the Planck scale, 
the photon energy ko has to appear as a second dimensionful scale. Then e.g. 

(Atf-i. (10.1) 



Because of fco = £ • k with £ = (1,0,0,0), a preferred spacetime direction comes into play 
which may lead to a violation of Lorentz invariance. 
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keep (1) (2a), (2b) (2c) (2d) (3) (4) Lorentz violation 

/ / / / / not needed absent 

/ / / / / / 

/ / / / / / 

/ / / / / (?) 

/ / / / not needed massive photon 



TABLE IV. The fate of Lorentz invariance where, in short, the assumptions taken are listed as follows. 
I: effective theory of Eq. (2.4), (2a)/(2b): homogeneous/isotropic defect distribution, (2c): dense defect 
distribution, (2d): infinite spacctimc volume, (3): modified photon momentum square <C mass of scalar field 
square, (4): momentum transfer from particles to defects suppressed by Planck scale. 



• ... an infinite spacetime volume: In this case the dispersion relation of the photon stays 
the same since gauge invariance is not violated. Note that in the simple model proposed 
no boundary conditions are set on any field and, therefore, the photon does not feel the 
finite volume. However, if the sprinkling procedure changes because of the finite volume the 
photon dispersion relation might also change. Due to calculational difficulties we have not 
been able to examine this case. 

• ... a dense distribution of defects: If the defects are assumed to be separated by a large 
distance, their distribution can no longer be considered to be continuous. The defect positions 
Xi will appear explicitly in the tensor structure of the photon field equation. This may lead 
to a modified dispersion relation of the photon that involves these positions. 

• ... an isotropic and homogeneous (random) defect distribution: If the distribution is dense 
but anisotropic it defines a preferred direction C in spacetime. Then £ also shows up in the 
tensor structure of the photon field equation. This may lead to a modified dispersion relation 
of the photon and, therefore, Lorentz violation. 

• ... the photon momentum squared <C mass of scalar particle squared: The photon gets a mass 
and electromagnetic waves are damped. We discard this case as it does not correspond to 
physical reality for photons with an energy much smaller than the Planck energy. Otherwise 
we would not be able to observe light from distant galaxies. 

It has not been possible to derive the modified photon dispersion relation for the more complicated 
cases mentioned above. This may be done in a future research project. Besides that, we could 
not rigorously demonstrate the influence of a finite spacetime volume as we did not impose any 
boundary conditions on the photon field. 

The final results are summarized in Table ITVl 
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A 


Wx/W2 


1 


0.701693 


0.1 


0.966576 


0.01 


0.996615 


0.004 


0.998645 



TABLE V. In the current table different values A of dimensionless defect distances are considered. For 



specific values the numbers of Eq. ( A.l I are computed and compared to each other. 



Appendix A: How good is the assumption of a dense defect distribution? 

In this section we consider four-dimensional Minkowski space with dimensionless coordinates 
t^. The numbers 



wi = > exp(h^l M ) , w 2 = Q exp(ii M l M ) , (A.l) 

shall be computed, where (1^) = (1, 1, 1, 1). The sum in w\ runs over n defects contained in a four- 
dimensional unit-hypercube EI whose edges are supposed to lie parallel to the axes of the coordinate 
system. The space diagonal of the hypercube runs from the point (0,0,0,0) to (1,1,1,1). The 
defects are assumed to lie at equal distances A. Hence, the hypercube contains (A -1 + l) 4 defects 
leading to the density g = (A -1 + l) 4 . The integral in W2 runs over the same hypercube, as 
well. We compare the values w\ to W2 for different defect spacings A. In principle w\ corresponds 
to W2 for infinitesimal defect separation. In Table [V] we see that the integral is already a good 
approximation for the dimensionless distance A = 0.01. 

We can now replace the dimensionless scalar product t^la by the scalar product of two di- 
mensionful quantities: the wave vector k 11 and the spatial four- vector x^. Then A corresponds to 
a product of a wave vector K and a distance AX in configuration space. Assuming the photon 
energy E = 1 TeV the dimensionless distance A is in accordance with the following AXj^] 

AX = — = — , 1/100 1Q , he « 2 • 10- 21 m . (A.2) 

K 10 12 eV • 1.602 • 10~ 19 J/eV v ; 

Thus the approximation we use is already very good even if the defect separation lies many orders 
of magnitude above the Planck length. 

Appendix B: Perturbative Feynman rules of the modified theory 



The following Feynman rules (B.la) - (B.le) directly follow from the action given by Eq. (2.1) 



k 



-i^A(k), A(k) = -^—, (B.la) 



* = _ iH(kh ffW = _i_)_ (B . lb) 



We assume the standard dispersion relation of the photon. 
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iX (0) e^k litl k 2 , e = -iX^e^ki^e . (B.lc) 



P 



X <_ =JdAGMp)GM- (B.ld) 



The first Feynman rule gives the photon propagator of the modified theory in Feynman gauge, 
which corresponds to the photon propagator of standard QED (in this gauge). The second gives 
the propagator of the scalar field <j>, where this has to be connected to a single defect. The third 
describes the interaction between cf> and the photon and the fourth the "defect vertex" for a finite 
number J\f of defects in a box. Performing the limit M i— > oo (consider Ass. (2a) - (2d)) combined 
with the integration over all Lorentz transformations to an upper cut-off leads to the fifth Feynman 
rule that we use to compute the one-loop correction in this limit 

^X^ =E 6 S^(p + q). (B.le) 

P Q 

So each eft line has to begin or end at one defect. If this is not the case one of the momenta has to 
be set to zero. 

Respecting Ass. (2a) - (2d) the following effective vertex can be introduced for the interaction 
of a photon with the defect via the </>-field. The latter has been integrated out at one-loop level 
resulting in 

v ^^0^v^ 7 = -cmr(k) , c = zb*x 2 Q, (B.if) 

k k 



with W^(k) given by Eq. (7.1b). 



Appendix C: Derivation of the Passarino Veltman reduction 



In this section the Passarino- Veltman reduction of the tensor integral I ga given by Eq. (5.3a) 



in Sec. VB will be presented in detail. 



V = / d d Q • (CI) 

J [q 2 -l/(6(°)) 2 +k] 2 [(k-q) 2 -fie] 

Performing the reduction we can neglect all ie. First, we contract I Qa with — 2k e and this leads to 
-2k%„= [d d q 



[q 2 - l/(bW) 2 ] Z (k-q) 2 

d d 1- f d d (k 2 +q 2 )qa (C 

r g2-l/(&(0))2] 2 J [ g 2_ 1 /( 6 (0))2] 2 ( A ._ g )2 • 
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Second, by contracting the latter result again with —2k a we obtain the quantity AK\ defined in 
Sec. EE 



4if x = Ak e k"I. 



gcr 



d d q 



ddq (k 2 + q 2 )(-2k 



[q 2 - l/(6(°)) 2 ] z (k-q) 2 
k- + <i 2 /•„, (k 2 +q 2 ) 2 



+ d' 



q 2 -l/(b(°)) 2 r (k-q) 



(C.3) 



Now we still have to decompose the integrals I3 and I4 into the master integrals of Eq. (5.8) 
k 2 + q 2 



h = d a q 



k 2 + 



[q 2 - l/(6(°)) 2 
' A f d d q 



k 2 / d d c 



1 



[q 2 - 



1 



(bW) 2 JJ [ga-l/(6W) 2 ] : 



+ d 



+ d 



1 



[g2 _ l/(6(0))5 



(C.4) 



A = / d 



A; 4 + 2&V + g 4 

[g2 _ l/(6(0))2] 2 (fc_ g )2 



fe 4 / d d Q 



+ / d a q 



[ q 2 -l/(b(0)) 2 Y(k-q)< 



l/(6(°)) 2 ] (Jfe-g) 5 



+ 2/t 2 / d d g 



l/(&(°)) 2 r(fc- g ) 2 



s/5 



i/ 6 



(C.5) 



d d q 



+ 



[g 2 - l/(6(0)) 2 ] (* - g) 2 (6W) 2 J H [ q 2_ i/(6(0))2] 2 (fc _ g) 



dV 



(C.6) 



d<V 



+- 



d d q 



(fc _ g )2 (6 (0) )2 y * [ g 2_ 1/(6 (0))2] 2 (A; _ g)1 



d d q 



(6(0)) 4 7 [g 2 -l/(6(0)) 2 ] 2 (A;-g)^ 



d d g 



+ 



(6(0))2 J " [ g 2 _ l/( 6 (0))2] (fc _ ^2 ( 6 (0))4 _/ * r ?2 _ 1/(6 (0))2] 2 (fc _ q) 



d d q 



(C.7) 



The contraction of I Qa with the metric tensor leads to K2 which was also defined in Sec. V B 



= / d d q 



[q 2 -l/{b(°)) 2 Y {k-qf 



d d q 



+ 



[ q 2 _ l/(6(0))2] (k - q) 2 (6(0)) 2 J *[ q 2_ l/( b (0))2] 2 (fc _ 9 ) 



dV 



(C.8) 



37 



Appendix D: Computation of scalar integrals 



Finally we evaluate the scalar integrals of Eq. (5.8) in Sec. VB We use dimensional regulariza- 
— 2e and later on we employ the reasonable redefinition 

7E + ln(4^). (D.l) 



tion with d = 4 — 2e and later on we employ the reasonable redefinition 
1 _ 1 

e e 

The simplest integral without any external momenta can be computed as follows 



An 



1 



-i(2^)"/dfi,/d 9?T ^ 



/(6(o))2 
u 

{0) 2 _ d T{l - d/2)r(rf/2) 
~ A ^> Y{d/2) {b > 2 



-i(2 / u) 4 -V- d / 2 (6(°)) 2 - d r ( 1 



i(2 M ) 2 ^ 2+? (6(°)) 2 ^ 1 )r(e-l) 



17T 



(6(0) ) 2 

The remaining integrals can be evaluated in the same way 



+ 0(e) 



Bo 0, 



1 



1 



(6(o) ) 2 ' (b(°)y 



i(2^/i) 4 - d / dft d / dg = = 

v w r(d/2) K ' 2r(2) 
i(2 / u) 4 -V- d / 2 (6(°)) 4 - d r (2 - 

i(2 / u) 2? 7r 2+? (6(°)) 2? r(e) 

1 



,(6(°))V 

5 K""W' ) =(2M4 "/ ddg (^T 



(2TT^- d J dx I d d q 



/(6(°)) 2 + ie) [(A; - g) 2 + ie] 
1 



(&(0))2 

{2irfi) 4 - d £ dx J d d q [{q - kx) 2 - M 2 ]~ 2 



(D.2) 



(D.3) 



+ [k 2 -2k.q+^ j)2 ]x + ie 



where 



-M 



k 2 x 2 +[k z + 



+ ie. 



This then leads to 



£ —k 



1 



(6(o))i 



= ivr" 



dx In 



(6(°)) 2 

fc 2 x 2 _ ^2 + l/( 6 (0))2) x + 1 /( 6 (0))2 _ ^ 



(D.4) 
(D.5) 

0( B ) 
(D.6) 
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As already mentioned, the Co-integral is both infrared and ultraviolet convergent. Hence we can 
set d = 4 at the beginning and we obtain 



C °( fc ' ' ( 6 (0))2' ' ( 6 (0))I 



dV 



1 



q 2 -l/(6(°)) 2 +ie) p-g) 2 +ie] 



r(2 + i) 
r(2)r(i) 



dxi / d^2 6(1 — xi — X2) 



x X\ 



2 / dx 2 



1 



(&(°)) 2 

1 — X 2 



+ ie I xi + [(k - q) 2 + ie] X2 



[{q - kx 2 ) 2 - M 2 ] 3 ' 



(D.7) 



with M 2 given by Eq. (|D.5l). Finally this results in 

1 



C -k,0, 



1 



0, 



(6(o)) 2 ' ' (6(0))2 



i-7r / dx 



1 — x 



k 2 x 2 - (k 2 + 1/(M°)) 2 ) x + l/(6(°)) 2 



Now we want to compute the remaining one-dimensional integrals 



dx In 



dx In 



k 2 x 2 _ ( k 2 + l/( 6 (0))2) + l/( 6 (0))2 _ i£ 



fc 2 x 2 _ (tf _ l/(6(0))2) 

x — ie 



/c 2 x 



ln(// ) + / dx ln(x) + / dx In 
io io 

ln(/i 2 ) - 1 + p [y\n(y) - y] ^2+1/(6(0) )2 _ ie 
ln(/i 2 ) - 1 



fc 2 - l/(6(°)) : 



1€ 



+ 



1 

k 2 



1 



In 



(6(0)) 2 V(^ (0) ) 



1 



ie 



(6C0)) 5 



In 



(&(°))V 



In ( -A-- + 
2 + 



(6(0))2 
P_l/(fo(0))2 



ie 



k 2 + 



k 2 



In 



(6(°)) 2 

fc2 + l/( & (0))2 

l/(6(°)) 2 



ie 



(D.8) 



(D.9) 



k 2 Ic 



dx 



1 — x 



an— >1— x 



A; 2 x 2 - (yfc 2 + l/(6(°)) 2 )x + l/(6( ))' 



ie 



dx 



1 



In 



z-l + l/[(6(°))fc] 
1 



if 



In ( x — 1 + 



(&(°))A;p 



ie 



[(&(°))A;p 



ie 



In 



[(&(°>)fc] : 



1 - ie 



In 



1 - (b^fk 2 



ie 



(D.10) 
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